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In the evaluation of weak interacting massive particles (WIMPs) detection rates, the WIMP- 
nucleus cross section is commonly described by using form factors extracted from charge distri- 
butions. In this work, we use different proton and neutron distributions taken from Hartree-Fock 
calculations. We study the effects of this choice on the total detection rates for six nuclei with 
different neutron excess, and taken from different regions of the nuclear chart. The use of different 
distributions for protons and neutrons becomes more important if isospin-dependent WIMP-nucleon 
interactions are considered. The need of distinct descriptions of proton and neutron densities reduces 
with the lowering of the detection energy thresholds. 
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The elastic collision with a nucleus is the main mechanism used by modern detectors to identify neutralinos, or 
more in general weakly interacting massive particles (WIMPs) [lH5| ■ The cross section describing this collision can be 
separated in a part related to the interaction between the WIMP and the single nucleon, and another one describing 
how the struck nucleon is inserted inside the nucleus @ . While the former term is one of the unknown which we aim 
to study in this type of investigations, the knowledge of the latter term is solid, and well grounded on many years of 
nuclear structure investigation. 
' The procedures commonly adopted substitute the nucleon distribution inside the nucleus with the charge distri- 
bution which is experimentally known with high spatial resolution because it is extracted from the large body of 
-H " accurate elastic electron scattering data. In specific applications the Helm form factor Q is usually adopted because 
of its simple analytic expression whose parameters have been chosen to provide a good description of the empirical 
charge distributions overall the nuclear chart, but for very light nuclei. An improvement with respect to this approach 
has been proposed in Ref. [§[ where it was suggested to consider directly the measured charge distributions from the 
compilation of Ref. Q . 

The use of the experimental charge distributions in the calculations of the WIMP-nucleus cross sections implies 
some simplifications of the problem. In the description of the cross section we consider the WIMP interacting with 
' the point-like nucleon, either proton or neutron. On the other hand, the nuclear charge density is mainly sensitive to 
I . the proton distributions folded with the electromagnetic proton form factor. 

In this article we propose a methodology to get over these simplifications. Our idea is to consider proton and 
neutron distributions generated by Hartree-Fock (HF) calculations. A work based on these ideas, and carried out 
■ within a relativistic-Hartree framework, has been proposed in Ref. [Io| . 

The effective nucleon-nucleon interactions used in HF calculations have been tuned to reproduce various ground 
state observables, among them also the charge root mean square (rms) radii of many isotopes. The remarkable success 
in describing elastic electron scattering data, and consequently charge distributions, indicates the reliability of these 
calculations in the description of the nuclear ground states. 

In the present work we use a specific implementation of the non-relativistic HF theory, based on a density-dependent 
and finite-range nucleon-nucleon effective interaction [Til ] . The performances of our approach have been tested against 
those of HF model which uses a zero-range interaction and those of a relativistic Hartree approach [12| . We found that 
the three types of calculations provide very similar descriptions of the ground state properties of spherical and closed 
shell nuclei. For this reason we arc confident that our results are not simply related to a specific implementation of a 
particular nuclear structure model, but they represent the general predictions of any microscopic mean-field theory. 

The basic ingredients of our model are presented in Sect. Ulland we show and discuss our results in Sect. Mil and 
in Sect. HVI we draw our conclusions. We give in the Appendix the Fourier-Bessel expansion coefficients of the proton 
and neutron distributions obtained in our calculations, for all the nuclei we have considered here. 
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II. THE MODEL 



We write the probability that, in a detector with Nt target nuclei, a nucleus, with A nucleons, recoils with an 
energy E R after an elastic collision with a WIMP as 

= N T £ a xA (v x )P(E R ,v x )$(v x )f(v x )dv x , (1) 

where v x is the velocity of the WIMP with respect to the detector and u m i n is the minimum velocity that produces 
a detectable recoil in the apparatus. In the above equation, we have indicated with <j x a the WIMP-nucleus elastic 
cross section, with P the probability that the nucleus after been struck by the WIMP with velocity v x recoils with 
energy En, with $ the flux of WIMPs with velocity v x , and with / the probability that the WIMP has velocity v x . 

We express the number of target nuclei as the product of the detector mass times the Avogadro's number divided 
by the target molecular weight, Md Na/At- The WIMPs flux is given by the WIMPs density times the velocity. We 
express the WIMPs density as 6 X /M x c 2 , where S x is the energy density of the WIMPs, M x is the WIMP mass and c 
indicates, as usual, the speed of the light. 

After the collision with a WIMP of velocity v x , the value of the nuclear recoil energy depends on the scattering 
angle. The nucleus may recoil with energies values between zero and 



where Ma is the mass of the target nucleus and [i x a = M x Ma/{M x + Ma) is the WIMP-target nucleus reduced 
mass. The energetics of the collision is such that the scattering is dominated by the s wave, therefore, because of the 
spherical symmetry of the scattered wave function, the probability P is independent on the scattering angle, and can 
be expressed as 

P(E R ,v x ) = ^4- ^ > (3) 

with < E R < E R , ma , x . 

For the evaluation of / we consider a Maxwell's distribution of the velocities centered at an average velocity vq. In 
this way, we obtain for the differential rate of recoil the expression 

dR(E R ) S x M d N A M A c 2 2 c 2 exph(f min /i; ) 2 ] 
dE R M x c 2 A T GxA 2ii\ a <* t-o^ U 

The total event rate, per unit of time, in a detector with energy detection threshold _E" t h is 

R(E th )= r^MdE R . (5) 
J Eth dE R 

The characteristics of the detector define Md and the mass of the recoiling nucleus. The unknown inputs of the 
calculations are the WIMP mass M x , the energy density S x , the average velocity vq, and the WIMP-nucleus cross 
section <j x a- 

The WIMP-nucleus cross section can be separated in spin-independent (SI) and spin-dependent (SD) terms the 
latter one active only in odd-even nuclei. In this work we are only interested in the SI term which we write as 

af A = t XP \ZF P (q) + g„ P N F n (q)\ 2 , (6) 

where Z and N — At — Z are, respectively, the number of protons and neutrons inside the target nucleus. In the above 
expression we have factorized the elastic WIMP-proton cross section a xp , and we have indicated with g np = a xn / a xp 
the ratio between the WIMP-neutron and proton cross sections. All the nuclear structure information is contained 
in the proton and neutron form factors F p (q) and F n (q) which depend on the modulus q = |q| = X /2M X E R of the 
momentum transferred by the WIMP to the nucleus. The form factors are defined as the Fourier transforms of the 
proton and neutron density distributions 

F a (q) = [ exp(iq-r)p Q (r)d 3 r = 4tt / j (qr) pjr) r 2 dr , (7) 



where a = p for protons and n for neutrons and jo indicates the zeroth-order spherical Bessel function. 
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In the above equation we have already considered that in our calculations the densities distributions have spherical 
symmetry. The expression ^ of the cross section, where the proton and neutron numbers are explicitly written, 
implies that the densities should be normalized as 

/>oo 

P P .n(r)d 3 r = 4tt / Pp . n (r)r 2 dr = 1 . (8) 
Jo 

In addition to these distributions we have considered the matter distribution 

1 

2 



Pm = T, (Pp + Pn) , (9) 



and the charge distribution, p c y L . The latter one is obtained by folding the proton distribution with the electromagnetic 
nucleon form factor. The densities p m and p c h have been used to calculate the form factors F m (q) and F c h{q) with 
expressions analogous to the ([7]). 

As we have already pointed out in the introduction, the procedure commonly adopted in the literature considers 
the form factors obtained by using the experimental charge density distributions. This implies the assumption that 
proton and neutron density distributions coincide with the charge distributions, therefore, in this case, instead of Eq. 
© the expression 

4a = °xpF 2 ch {<D \Z + 9n P N\ 2 , (10) 

is used. 

In the present work, the proton and neutron distributions have been calculated by using a mean-field approach 
where the densities are given by the expression 

P(r) = E l&(r)| a ■ (H) 

k<k F 

In the above expression (f>k indicates the single particle (s.p.) wave function characterized by the quantum numbers 
identified by k and the sum is limited to all the states below the Fermi surface. Clearly, proton and neutron densities 
are obtained by selecting in the sum only those s.p. states related to protons or neutrons respectively. 
In our calculations the s.p. wave functions are obtained by solving the HF equations [13j 

- ^ V 2 ^(r) + EWfc(r) - J dVW(r,r')<fe(r') = e fe <fe(r) , (12) 
where m is the nucleon mass, the s.p. energy and we have indicated the Hartree term as 

U(r) = J dW(r,r>(r') , (13) 

and the Fock - Dirac term as 

W(v,v')= £ y(r,r')^(r')^(r) ■ (14) 

k' <k F 

In the above equations, V(r, r') represents the effective nucleon-nucleon interaction. In our calculations we use an 
interaction which has finite-range character in the scalar, isospin, spin and spin-isospin terms and zero-range character 
in the spin-orbit and density-dependent terms. This interaction, known in the literature as Gogny interaction [llj . 
contains 13 parameters whose values have been chosen to reproduce experimental masses and charge radii of a large 
number of nuclei. We have done our calculations with two different choices of the parameter values, the more 
traditional D1S force [3], and the more modern DIM force [HI]. Since we did not find relevant differences in the 
results obtained with the two effective interactions, we show here only the results obtained with the DIM force. 

In our calculations we treat spherical nuclei, therefore we found convenient to write the s.p. wave functions by 
factorizing a term R n ij(r) depending on the distance r from the nuclear center and another part related to the 
spherical harmonics. In this case, the density (fTTj) assumes the expression 



= 7-E u ny(2i + l) |iWr)| 2 ■ (15) 
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BE/nucleon [MeV] rms [fm] 



Nucleus 


HF 


Exp. 


HF 


BCS Exp. 


ib 


7.98 


7.98 


2.76 


- 2.74 


40 Ar 


9.32 


8.60 


3.40 


3.40 3.42 


40 Ca 


8.51 


8.55 


3.47 


- 3.48 


72 Ge 


8.57 


8.73 


4.02 


4.02 4.05 


136 Xe 


8.44 


8.40 


4.80 


4.81 - 


208p b 


7.83 


7.87 


5.48 


- 5.50 



Table I: Binding energies (BE) per nucleon and charge rms radii calculated with the DIM interaction by using the HF and 
BCS models and compared with their experimental values. 



where n, I and j indicate, respectively, the principal quantum number, the orbital and total angular momentum 
characterizing the s.p. state whose degeneracy is 2j + 1. In Eq. (fPo)) we have indicated with w^ . the occupation 
probability of the s.p. state. In closed shell nuclei we have v 2 y = 1 for the states below the Fermi surface, and 
v nij = f° r those above it. In our calculations we have studied three nuclei of this type, 16 0, 40 Ca and 208 Pb which 
are representative of various regions of the nuclear chart. 

The performances of our calculations in reproducing ground state observables are rather good as it is shown in 
Table|T]where we compare binding energies and charge rms radii with their experimental values. The good description 
of the ground state properties within the HF theory is not related to the present implementation of the mean-field 
approach, but it is a feature of this type of calculations, as it is shown in Ref. [l2j], where our results are compared 
with those obtained by using a zero-range interaction and also with those produced within a relativistic framework. 

Because of their use in WIMPS detectors HU, we have also considered three open shell spherical nuclei, 40 Ar, 
72 Ge and 136 Xe. In the HF picture, in these nuclei the s.p. levels are not fully occupied. In this case, the pairing 
phenomena, irrelevant in closed-shell nuclei, cannot be neglected. We have considered pairing effects by solving the 
equations of the Bardeen, Cooper and Schriefer (BCS) theory applied to finite nuclear systems [13j . We start from 
the s.p. wave functions and energies generated by our HF calculations, and then we solve the set of equations: 

A a = -— 1— V y/2j b + lu b v b (bW\V\aaO) , (16) 
2^/2j a + 1 ^ 




In the above equations we have used u 2 = 1 — v 2 , a and b indicate different sets of the n,l,j quantum numbers, 
and N is the total number of nucleons subject to the action of the pairing. The matrix element of the interaction 
is calculated by coupling to zero the angular momenta of the a and b s.p. states. The solution of these equations 
provides the values of A which is the gap energy, i. e. the energy difference between the last occupied level and the 
first empty level in the HF picture. The other output of the calculation, i. e. the values of the v a , is more relevant 
for our purpose since, as expressed by Eq. (|15p . it modifies the density distributions. 

We carried out the BCS calculations by using the same nucleon-nucleon interaction used in the HF calculations, 
but without Coulomb and spin-orbit terms, as it is commonly done in the literature for this type of calculations. The 
relevance of the pairing can be deduced by observing how much the values of v differ from those of the HF picture. 
We show in Table ITll the values of v 2 related to the partially occupied s.p. states nearby the Fermi level, for the three 
open shell nuclei we have considered. The values of v 2 for the neutrons in the 136 Xe nucleus clearly indicate the shell 
closure for N=82. 



III. RESULTS 



We have calculated proton, neutron and charge density distributions within the theoretical framework outlined in 
the previous section. The charge distributions have been obtained from the proton distributions by folding them 
with the proton electric form factor extracted from elastic electron-proton data. The results we show here have been 
obtained by usin g a dipole parametrization of the proton form factor [l7j ] . We have verified that more accurate 
parametrizations [18j modify our charge distributions on few parts on a thousand. 
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v 2 




Nucleus 


protons 


neutrons 


4U Ar 


1^3/2 I/7/2 
u.uou u.uuo 


2P3/2 

0.003 


I/7/2 
0.257 


r2 Ge 


2P3/2 I/5/2 
0.541 0.286 


2pi/ 2 
0.984 


I39/2 
0.013 


m Xe 


l<?7/2 2d 5 / 2 

0.388 0.115 


1.000 


lftg/2 
0.000 



Table II: Values of v 2 for the partially occupied s.p. levels near the Fermi level. 



In Fig. Q]we compare the charge distributions obtained in our HF and HF+BCS calculations with the empirical 
ones extracted from elastic electron scattering experiments. In this figure we present the distributions multiplied 
by r 2 to directly show the functions which are integrated in Eq. ([7]). We observe a very good agreement of our 
results with the experimental distributions. This is a consequence of the fact that the charge rms radii are some of 
the experimental data included in the fit procedure used to define the parameters of the DIM force. As it is well 
known, the main differences between experimental and theoretical charge distributions arise in the nuclear interior. 
In general, the experimental charge distributions are smoother than the theoretical ones. A well known example of 
these phenomena is the case of 208 Pb which has been widely investigated from both experimental [19( and theoretical 
[20I l21j points of view. 

For the open shell nuclei, we compare the results obtained in both HF and HF+BCS frameworks, dashed and dotted 
lines, respectively. For the two nuclei where the empirical charge distributions are available, nuclei the HF+BCS 
distributions provide a better description of the experimental ones than those obtained in the HF calculation. This 
is more evident in the case of the 72 Ge nucleus. 

The results of Fig. Q] provide an indication of the capacity of our model to describe the charge distributions. 
The WIMP-nucleus interaction is active on the full matter distribution which we claim to describe with the same 
accuracy than the charge distribution. Since no accurate data on matter distributions are available we cannot test 
this assumption, which, on the other hand, we believe it is well grounded. 

In the literature, properly normalized charge distributions are used to estimate WIMP-nucleus cross sections instead 




Figure 1: Charge distributions obtained with our HF (dashed lines) and HF+BCS (dotted lines) calculations compared with 
the empirical distributions (full lines) [3], extracted from elastic electron scattering data. The distributions, normalized as 
indicated by Eq. (|8]), have been multiplied by r 2 . 
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Figure 2: Differences between proton and charge distributions (dotted lines) and between matter and charge distributions 
(dashed lines), Eq. 1|18[) . multiplied by r 2 . 



than matter distribution. In order to estimate the error done in this approximation we have calculated the difference 
between proton (a = p) and matter (a — m) distributions with respect to the charge distributions 

Ap a {r) = p a (r) - p ch {r) . (18) 

We show in Fig. [5] these differences, multiplied by r 2 , obtained in the HF framework for the 16 0, 40 Ca and 208 Pb 
nuclei, and in the HF+BCS framework for the 40 Ar, 72 Ge and 136 Xe nuclei. Unless stated otherwise, the results 
presented henceforth have been obtained in this manner. 

The first remark about the results of Fig. [2] is that Ap p {r) and Ap m (r) have the same order of magnitude in all 
the nuclei considered. This indicates that the contribution of the electromagnetic proton form factor is not negligible. 
In the opposite case, we would have found smaller differences with the proton distributions than with the matter 
distributions. A detailed discussion on the WIMP-nucleon interaction is not in the scope of the present article. We 
point out, however, that in the expression ^ of the cross section, the WIMP-proton cross section is factorized. From 
the theoretical point of view, this implies the requirement of using point-like densities to calculate the form factors 
F p (q) and F n (q) in order to avoid a kind of double counting between electromagnetic and weak nucleon form factors. 

A second remark concerning the results shown in Fig. [5] indicates that in nuclei with equal number of protons 
and neutrons, 16 and 40 Ca, the proton and matter distributions, when normalized as indicated in Eq. ([5]), are very 
similar. The differences start to appear in the 40 Ar nucleus where Z — 18 and N = 22. The larger is the difference 
between Z and N, the larger is the difference between proton and matter distributions. In 72 Ge, 136 Xe and 208 Pb 
nuclei Ap p and Ap m differ mainly in the nuclear interior. In this region, the differences with the charge distributions, 
here considered as reference result, are even out of phase. The situation becomes more stable in the surface region 
where the curves show similar behaviors. The negative values of all the curves in the surface region is due to the 
fact that the charge distributions are always wider in space than the point-like distributions. This is the effect of the 
folding with the proton electromagnetic form factor. Interesting to notice that in all the cases, with the exceptions 
of 16 and 40 Ca, the Ap p is always larger than Ap m . Our calculations indicates that the identification of the charge 
distributions with the full matter distributions is a better approximation than that obtained by identifying charge 
with proton distributions. 
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Figure 3: (Color on line) Nuclear form factors calculated with various density distributions. Full, dotted and dashed lines show 
the results obtained by using, respectively, the charge, proton and matter distributions obtained in our calculations. The red 
dashed-doubly-dotted lines are the form factors obtained by using the empirical charge distributions. The dashed-dotted lines 
indicate the form factors obtained within the Helm model. The vertical blue lines indicate the value of the momentum transfer 
corresponding to a nucleus kinetic energy of 100 keV. 



The form factors to be used in Eq. ((6]) to calculate the WIMP-nucleus cross section are shown in Fig. [3] where we 
have indicated, respectively, with full, dotted and dashed lines the results obtained by using the charge, proton and 
matter distributions obtained in our calculations. In the figure we also show with the red dashed- doubly-dotted lines 
the form factors obtained by using the empirical charge distributions. The dashed-dotted lines show the Helm form 
factors [7j] which are commonly used in the literature to calculate the WIMP-nucleus cross section. The Helm form 
factor is obtained by considering a charge density with a Gaussian surface distribution, and has a simple analytic 
expression 

\F H clnM)\ 2 = 

where j\ is the first-order spherical Bessel function, and R and s are two parameters which are chosen to reproduce 
at best the experimental charge density distributions on all the nuclear chart. We used the parametrization suggested 
in Refs. 0, ill 

R 2 = c 2 - bs 2 + 18.65 fm 2 , (20) 
c = fl.23 A 1 / 3 - 0.6) fm , (21) 
s = 0.9 fm . (22) 

The differences between the various form factors increase with the value of q. This reflects the fact that, by 
increasing the resolution power, the probe is able to perceive the differences between the various densities. The form 
factors drop rather quickly with increasing value of q. The relevance of these differences in the calculation of the 
total event rate can be estimated by considering that, in the expression ([5|), the integration starts from the value of 
the threshold detection energy E t h- As a reference example, in Fig. [3] we have indicated with vertical blue lines the 
values of the momentum transfer corresponding to a recoil energy, Er, of 100 kcV. This is a large value with respect 
to the performances of modern detectors. In this example, the form factors to the left hand side of the vertical lines 
do would not contribute to the total rate. This indicates that the largest values of the form factors are excluded, 
confirming an obvious consideration, the lower is the detection energy threshold, the higher is the rate of the detected 
events. More interesting is the comparison between the various nuclei. The heavier is the nucleus, the higher is the 



qR 



cxp(-gV) , 



(19) 
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Figure 4: Differences 5 a (Eu), as defined in Eq. (|23p , between the differential rates {TJ calculated with our density distributions 
and those obtained with the Helm form factors. The values for the charge (solid lines), proton (dotted lines) and matter (dashed 
lines) distributions are shown. In each panel, the multiplication factors used to enlarge the y scale are indicated. 



value of q required to obtain a 100 keV recoil energy. For the two heavier nuclei 136 Xe and 208 Pb, the lines appear 
after the first diffraction minimum, therefore, for these nuclei, the differences between the various densities can be 
more relevant. 

We used our density distributions to calculate the differential, and total event rates, Eqs. |T]) and @. The values 
of the input parameters have been chosen coherently with what is used in the literature. If not stated otherwise, the 
results of the calculations we present in this paper have been obtained by assuming a xp = 10~ 42 cm 2 = 10 -16 fm 2 , 
S x = 0.3 GeV/fm 3 , v = 0.001 c, g np = 1, M x = 100 GeV. 

The first results we want to discuss are related to the differential rates of eq. (H|). We show in Fig. 0]the differences 



~dR(E R y 




~dR(E R y 




dE R 


a 


[ dE R 


Holm 



(23) 



between the differential rates obtained with the distributions found in our calculations and those obtained by using 
the Helm form factors as a function of the recoil energy E R . In these calculations we used a threshold detection 
energy of 10 eV only. Full, dotted and dashed lines show the differences obtained by using the charge (a = ch), 
proton (a = p) and matter (a = m) densities, respectively. 

A first observation is that these differences are rather small. The relative differences in the peaks of the figure 
are 5% at most. A second observation is that these differences increase with increasing mass number. The order of 
magnitude of these differences in 208 Pb is 1000 times larger than that of 16 0, 100 times larger than those of 40 Ca, 
40 Ar and 72 Ge and, only, 10 times larger than that of 136 Xe. In all the nuclei we have considered, with the exception 
of 40 Ca, the smallest differences we found are those with the charge distributions. This could be expected since the 
parameters of the Helm form factors have been chosen to describe charge distributions. For the nuclei with Z ^ N 
the results obtained with the proton distributions differ more than those obtained with the matter distributions. This 
is a direct consequence of what we have observed in Fig. [2j 
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Figure 5: Total rates, as given by Eq. as a function of the threshold detection energy. The values found for the charge 
(solid lines), proton (dotted lines) and matter (dashed lines) distributions as well as those obtained with the Helm form factors 
(dashed-dotted lines) are shown. The values of the 16 O nucleus have been multiplied by 10. 



We show in Fig. [5] the total rates R a (E t h) given by Eq. (|5|) as a function of the threshold energy. The different 
lines indicate the results obtained with the various densities. The behavior of the various lines is well understood. By 
increasing the value of the threshold energy, which means a lowering of the sensitivity of the detector, the number of 
the events becomes smaller. The comparison of the results obtained for the various nuclei indicates that at low values 
of E t h the heavier nuclei are more efficients in detecting events, because of the larger number of target nucleons inside 
the nucleus. We point out that, in the figure, the 16 results have been multiplied by a factor 10. On the other hand, 
the total rates drop more quickly in heavier nuclei than in lighter ones. At 40 keV the total rates expected in 208 Pb 
are much smaller than those expected in the medium-heavy nuclei we have considered, 40 Ca, 40 Ar and 72 Ge. 

The scale of Fig. [5] does not allow to appreciate the differences between the various results. For this reason, we have 
calculated the relative differences between the total rates obtained by using the matter distributions with respect to 
the results obtained with the Helm form factors 

An fzp \ RrajEth) - ggelmCgth) (OA s 

/\H m (E th ) = — — . (24) 

^Helmi-C'thJ 

These relative differences are shown in Fig. [5] as a function of E t h- In this figure, in addition to the results related to 
total rates presented in Fig. [5j we show also results obtained by changing the WIMP mass M x . 

The value of the relative differences increases the heavier is the nucleus. We obtain a maximum value of 3% in 16 
and 10-20% in 208 Pb. The general enhancement with increasing Eth is more related to the lowering of the value in 
the denominator of Eq. (|2~4"|) than to a real increase of the difference between the results obtained with matter and 
Helm distributions. 

By comparing the results obtained with different values of M x , we observe that AR increases with increasing mass. 
Since the minimum value of the momentum transfer for the detection is q m in — v/2£th M x , if M x increases, also g m ; n 
increases. Therefore, a larger part of the form factors at low q, shown in Fig. [31 is excluded from the integral of Eq. 
([2]), and the final result is more sensitive to the differences between various calculations, differences which appear at 
high q values. In heavier nuclei, the increase of E t h moves the limit of this integration nearby the first minimum of the 
form factor, which is slightly displaced in the Helm model with respect to that obtained with the matter distribution. 

To estimate the relevance of the differences in the use of the various form factors, we have calculated hypothetical 
values of the WIMP-proton cross section for which a rate of a single detected event per day is obtained in a 100 ton 
detector. For all the nuclei we have investigated, we show in Fig. [7] these limits, obtained for Eth = 1, 10 and 100 
keV, as a function of the WIMP mass. The dashed-dotted curves correspond to the values found with the Helm form 
factors, while the solid lines indicate, respectively, those obtained by using our HF (for 16 0, 40 Ca and 208 Pb) and 
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Figure 7: Upper limits of the WIMP-proton cross section, a xp , as a function of M x , for the detection of a single event per day 
identified in an ideal detector of 100 tons, for Eth = 1, 10 and 100 keV. The solid lines show the results obtained by using our 
matter distributions obtained in HF or HF+BCS calculations, while the dashed-dotted lines indicate the results obtained by 
using the Helm form factors. The E t h = 100 keV results have been multiplied by 0.1 in case of ls O, 40 Ca, 40 Ar, 72 Ge nuclei, 
and by 0.01 for 136 Xe and 208 Pb nuclei. 
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HF+BCS (for 40 Ar, 72 Ge and 136 Xe) matter distributions and HF+BCS (for 40 Ar, 72 Ge and 136 Xe). The E th = 100 
keV results have been multiplied by the factors indicated. 

Obviously, the lower is Eth, the lower is the line of the exclusion plot. The differences between our results and those 
obtained with the Helm form factors become larger with increasing target mass, with increasing M x values, and with 
increasing Eth- These results show a minimum difference of 1.4% in the case of 16 for Eth = 1 keV, and a maximum 
one of 15% for 208 Pb with Eth = 100 keV. 

The results we have so far presented have been obtained by assuming the same coupling strength of the WIMP 
with protons and neutrons, i. e. g np — 1 in Eq. ([5]). We expect that, by releasing this assumption, as suggested, 
for example, in Ref. [23| . the requirement of a correct description of proton and neutron distributions becomes more 
important. For this reason, we have carried on calculations for different values of g np . We have chosen values of g np 
producing extreme situations. The largest and smallest values we have considered, ±2, slightly exceed those indicated 
in [HI, values based on a compatibility analysis of the DAMA Q, H], CoGeNT and XENON100 0] data. The 
other values of g np have been chosen because of their physical meaning. The isospin conserving coupling, which has 
been our reference up to now, is obtained with g np = 1. The value g np = selects only the WIMP-protons interaction. 
Finally, by using g np = — 1 we generate a cross section sensitive only to the differences between proton and neutron 
distributions. 
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Figure 8: Total rates, Eq. ([5]), calculated for different values of g np , which are indicated in the insert of the 16 panel. 



The total rates calculated with these values of g np , by using M x = 100 GeV, are shown in Fig. [8] In all the 
calculations the sequence of the various responses is similar. The lower lines are obtained with g np = — 1 the situation 
where the proton and neutron contribution cancel with each other. In the calculations where F p = F n , for example 
those which use the Helm form factors, the total rates for 16 and 40 Ca are exactly zero. Next we have the results 
with g np = where the WIMP interacts only with the proton. We have then the results with g np = —2 where the 
WIMP-neutron interaction dominates. In nuclei with N = Z the total rates almost ove rlap with that obtained with 
gnp = 0. For positive values of g np the total rates increase with increasing g np . In Ref. [231 ] a value of about -1.5 has 
been suggested. 

The aim of our work is to study the need of using accurate proton and neutron distributions to describe the WIMP- 
nucleus cross section. For this reason, we show in Fig. |9]the relative differences AR, Eq. (f24|) . obtained with the 
various values of g np . In the two N — Z nuclei, 16 and 40 Ca, the results with g np = — 1 are not indicated since in 
this case i?Hcim is zero with the obvious consequences in the calculation of AR. 

The behavior of the various lines has different characteristics for the N = Z nuclei and for the other ones. In the 
former cases the smaller differences appear for g np = and the largest ones for g np = — 2. The relative differences 
can reach values around the 20% in heavy nuclei. 

We have evaluated the consequences of these differences by calculating limits as those shown in Fig. [7] with 
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Figure 9: Relative differences between total rates cafculated with our matter distributions and those obtained with the Helm 
form factor for M x — 100 GeV and for different values of g np , indicated in the insert of the upper left panel. 



different values of g np . We show in Fig. [TU] the relative differences Aer xp between the limits of the WIMP-proton 
cross section calculated with our matter distributions and those obtained with the Helm form factors. The results for 
E t h = 1) 10, 100 keV are shown. The black (full, dashed-dotted and short-dashed) lines show the relative differences 
of our reference calculations done with g np = 1, i.e. those obtained from the results of Fig. [JJ The red (dotted, 
dashed-doubly-dotted and log-dashed) lines indicate the differences obtained with g np = — 2, a value comparable with 
that suggested in Ref. [23]. In all the cases the largest differences have been obtained for a threshold energy of 
100 keV. As it has been previously discussed, the differences arc enhanced with increasing Eth.- We do not identify 
common trends in this figure. Sometime the g np — —2 results show the largest differences, in other cases the largest 
differences are produced by g np = 1. We observe that in 40 Ar, 72 Ge and 136 Xe these values can reach the 20%. 



IV. CONCLUSIONS 



In this work, we have studied the validity of the approximation commonly done in the literature consisting in 
considering the nuclear charge density instead than the matter density in the calculation of the WIMP-nucleus cross 
sections. While charge distributions are experimentally known for various nuclei with high accuracy Q, the matter 
distributions are essentially unknown. For this reason, we propose to use the point-like matter distributions obtained 
in mean-field calculations of HF type. We present here the results obtained with a specific implementation of the 
HF calculations, involving an effective finite range interaction of Gogny type [ll[. Our experience [l2[ indicates that 
these results are general and more related to the mean-field description of the nucleus rather than to its specific 
implementation. 

We have shown that our model describes well the experimental charge distributions, and we have assumed that it 
is also able to properly describe the matter distributions. We have compared the form factors calculated with our 
matter distributions with the Helm form factors commonly adopted in the literature. The differences between the 
factors calculated by using different distributions show up at high values of the momentum transfer, usually after the 
first diffraction minimum. The relevance of these differences is strictly related to the value of the detection threshold 
energy E^, and decreases with it. 

We have calculated differential, total rates and upper detections limits. In this last case, we found that the differences 
between our results and those obtained with the Helm form factor become larger with increasing target mass, with 
increasing M x values, and with increasing E t h- We found a maximum relative difference of about 15% for 208 Pb with 
E th = 100 keV. 

The requirement of using matter instead than charge distributions in the calculation of the WIMP-nucleus cross 



13 




-40 ' 1 ' 1 ' 1 1 1 ' 1 1 -20 

100 200 300 400 500 



100 200 300 400 500 



M % [GeV] 



M x [GeV] 



Figure 10: (Color on line) Relative differences between the limit values of a xp calculated by using our matter distributions and 
those obtained with the Helm form factors. The results for E t h = 1 (full and red dotted lines), 10 (dashed-dotted and red 
dashed-double-dotted lines) and 100 keV (short-dashed and red long-dashed lines) are shown. Black full, dashed-dotted and 
short-dashed lines correspond to g np — 1.0, while red dotted, dashed-double-dotted and long-dashed lines correspond to the 
isotope violating hypothesis with g np — —2.0. 



section becomes more important when the assumption of isospin conserving WIMP-nucleon interaction is released. 
In this case, the proton and neutron distributions are mixed with different weights, therefore an accurate description 
of the two distributions is mandatory. Also in this case, the value of E t h is the parameter which dominates the value 
of the uncertainty, however, the use of mixing values close to those suggested in the literature increases the difference 
with respect the results obtained with the Helm model. In specific situations the relative differences can reach the 
value of 20%. 

In the estimation of the total event rate related to the WIMPs detection, the values of many input quantities are 
unknown, and strong assumptions are done on them. On the other hand, the nuclear physics part of the process is 
well understood, and we propose to use the modern nuclear structure results to improve its description. 



Appendix: Fourier-Bessel coefficients of density distributions 



In this appendix we show the values of the Fourier-Bessel coefficients which allow the reconstruction of our proton 
and neutron density distributions for the nuclei we have considered. These distributions can be obtained by the 
expression 

^r) = E^, (25) 

where the are defined as 

with R the maximum value of the radius where the density is supposed to be different from zero. The density 
distributions obtained with the coefficients given in Tables HTT1 and HVt have the usual nuclear physics normalizations, 
i. e. the integrals of Eq. {[8J) are normalized to the number of protons or neutrons. 
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■o 



'Ca 



proton 



neutron 



proton 



neutron 



proton 



neutron 



1 2.89514E- 

2 5.44734E- 

3 2.37463E- 

4 -1.51728E 

5 -1.94164E- 

6 -5.40039E 

7 1.85720E- 

8 6.24424E- 

9 -1.18242E- 

10 -9.50609E 

11 -1.71403E 

12 4.61957E- 

13 -1.93572E 

14 -7.29176E 
15 



02 2.90605E-02 
02 5.53342E-02 

02 2.50485E-02 
-02 -1.46814E-02 
-02 -1.99175E-02 
-03 -6.14889E-03 

03 1.34614E-03 

04 3.75985E-04 
-03 -1.29026E-03 
■04 -9.80391E-04 
-04 -1.66395E-04 

05 6.54090E-05 
-05 -9.66199E-06 
-05 -6.72835E-05 



6.17474E- 
5.75764E- 
-2.94985E 
-1.93097E 
1.87819E- 
9.02298E- 
-1.89926E 
1.40503E- 
2.91761E- 
1.34965E- 
1.77650E- 
1.68542E- 
1.13150E- 
2.12233E- 
2.61156E- 



02 6.23662E-02 
02 6.08146E-02 
-02 -2.82238E-02 
-02 -2.10479E-02 

02 1.85152E-02 

03 9.93418E-03 
-03 -2.77723E-03 
03 -1.48915E-03 
03 -1.02014E-03 
03 -2.77186E-03 
03 -2.38989E-03 
03 -1.69061E-03 

03 -8.69682E-04 

04 -4.41138E-04 
04 4.04638E-05 



7.69622E- 
1.87236E- 

-6.54988E- 
3.20906E- 
1.42479E- 

-2.81145E- 
1.43845E- 
1.53527E- 
3.00274E- 
3.76239E- 

-4.27147E- 

-3.33174E- 
3.06557E- 
1.80972E- 

-3.10532E- 



02 1.16349E-01 
02 2.73910E-02 
-02 -8.28045E-02 
02 4.81394E-02 
02 1.17763E-02 
-02 -2.99372E-02 
02 1.88786E-02 

02 -6.90422E-03 

03 -2.34458E-02 
03 -1.08603E-03 
-03 5.77452E-03 
-03 2.14901E-03 
03 -8.57191E-04 
03 -3.12159E-03 
■04 4.46753E-04 



R [fmj 


7.0 


7.0 


7.0 


7.0 


10.0 


10.0 



Table III: Fourier-Bessel coefficients A M and radii R, see Eq. (|25[) , f° r the proton and neutron densities obtained in the HF 
calculations. 







4u Ar 


72 Ge 


136 Xc 




proton neutron 


proton neutron 


proton neutron 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



3.09708E- 
5.88325E- 
2.10665E- 
-2.71647E 
-2.58345E 
1.35027E- 
1.14853E- 
3.83701E- 
-2.49991E 
-2.09091E 
-1.66806E 
3.19675E- 
8.86515E- 



02 3.76246E-02 
02 7.02530E-02 

02 2.43967E-02 
-02 -3.07229E-02 
-02 -2.67849E-02 

03 5.39599E-03 

02 1.73114E-02 

03 9.32738E-03 
-03 2.11831E-03 
■03 1.01436E-03 
-04 1.16466E-03 

04 4.98071E-04 

05 -8.27115E-05 



4.98732E 
6.19545E- 
-1.96985E 
-3.69138E 
1.29234E- 
2.33454E- 
3.16112E- 
-8.94086E- 
3.68467E- 
2.56986E- 
1.12422E- 
-2.13356E 
7.99035E- 
7.52024E- 
4.16595E- 



02 
02 
-02 
-02 
02 
02 
03 
-04 
03 
03 
04 
-04 
05 
05 
05 



18996E-02 
60695E-02 
.88170E-02 
.41701E-02 
68135E-02 
51003E-02 
.37351E-02 
.58768E-02 
.21935E-03 
.03807E-03 
28180E-04 
.48872E-04 
73510E-06 
45816E-05 
51430E-05 



4.64213E- 
5.66577E- 
-2.89511E 
-4.44567E 
2.16418E- 
2.99702E- 
-1.22571E 
-2.07746E 
-4.58139E 
-8.95531E 
-3.60592E 
-1.72562E 
4.85598E- 
4.04238E- 



02 6.96034E-02 
02 8.15343E-02 
-02 -3.93097E-02 
-02 -5.02359E-02 
02 3.33255E-02 
02 2.78136E-02 
-02 -2.05113E-02 
-02 -7.97235E-03 
-03 2.24455E-02 
-04 1.91808E-02 
-03 5.18993E-03 
-03 1.49995E-03 
04 1.63508E-03 
04 6.51292E-04 



R [fmj 


9.0 


9.0 


9.0 


9.0 


11.0 


11.0 



Table IV: Fourier-Bessel coefficients and radii R, see Eq. (|25[) , for the proton and neutron densities obtained in the HF + 
BCS calculations. 
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